Thermodynamical description of stationary, asymptotically 
flat solutions with conical singularities 



Carlos Herdeiro/'El Eugen Radu, 2, and Carmen Rebelo 1 '!! 

1 Centra de Fisica do Porto - CFP 
Departamento de Fisica e Astronomia 
Faculdade de Ciincias da Universidade do Porto - FCUP 
Rua do Campo Alegre, ^169-007 Porto, Portugal 
2 Institut fur Physik, Universitdt Oldenburg 
Postfach 2503 D-26111 Oldenburg, Germany 
(Dated: April 2010) 

We examine the thermodynamical properties of a number of asymptotically flat, stationary (but 
not static) solutions having conical singularities, with both connected and non-connected event hori- 
zons, using the thermodynamical description recently proposed in [J. The examples considered are 
the double-Kerr solution, the black ring rotating in either S 2 or S 1 and the black Saturn, where the 
balance condition is not imposed for the latter two solutions. We show that not only the Bekenstein- 
Hawking area law is recovered from the thermodynamical description but also the thermodynamical 
angular momentum is the ADM angular momentum. We also analyse the thermodynamical stability 
and show that, for all these solutions, either the isothermal moment of inertia or the specific heat 
at constant angular momentum is negative, at any point in parameter space. Therefore, all these 
solutions are thermodynamically unstable in the grand canonical ensemble. 

PACS numbers: 04.20.-q, 04.20.Dw, 04.20.Jb, 04.50.Gh, 04.70.Dy 



1. INTRODUCTION 

The existence of conical singularities on a manifold implies geodesic incompleteness and there- 
fore these singularities must be regarded as space-time boundaries from the viewpoint of classical 
differential geometry. Semi-classically, however, solutions with conical singularities havegeneri- 
cally a well defined Euclidean action, and thus well defined thermodynamical properties |2j ■ This 
suggests that such solutions may be associated to well defined quantum states in an appropriate 
theory of quantum gravity, upon a correct quantisation of the physical parameters, including the 
conical deficit/excess. Indeed there are special examples of manifolds with conical singularities 
in which a quantum theory of gravity is well defined. This is the case of string theory on certain 
orbifold d, [j|; m these backgrounds string theory is solvable and smooth, due to new degrees of 
freedom localised at the tip of the cone - the twisted sector. 

Even at the classical level, there is usually a clear physical interpretation for the conical 
singularity. For instance, as providing an otherwise impossible mechanical balance between two 
black holes. In this class of examples, studying the singularity's properties unveils features 
about black hole interactions, as shown in 0, Q in the case of the double Kerr solution in 
four spacetime dimensions. Thus, despite their incomplete nature, backgrounds with conical 
singularities have well defined physical properties worth investigating. 
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In this paper we shall continue to investigate the thermodynamical description introduced in 
[l[ . The basic observation therein was that a natural choice of thermodynamical variables leads 
to the Bekenstein-Hawking area law for the entropy, even for solutions with non-connected event 
horizons, in contrast to previous approaches. All the examples given in [l| were static solutions. 
Here we consider various stationary solutions which, similarly to are asymptotically flat, i.e 
the conical singularity does not extend to spatial infinity. We show that not only the area law 
is obtained for the entropy, but also that the angular momentum computed thermodynamically 
from the free energy coincides with the ADM angular momentum, unlike the thermodynamical 
description in [||. Since we have a consistent thermodynamical description of these solutions, we 
can also examine their thermodynamical stability. The generic feature found is that they are all 
unstable in the grand canonical ensemble. 

This paper is organised as follows. In section 2 we introduce the appropriate thermodynamical 
variables, the various energy functionals and all relevant formulas to be applied to the specific 
examples. These are considered in the following sections: black rings in section 3, the double 
Kerr in section 4 and the black Saturn in section 5. We present some final remarks in section 6. 

2. FORMALISM 

In order to discuss the thermodynamical properties of such backgrounds an appropriate set 
of thermodynamical variables must be chosen. It was suggested in [l[ that the first law of 
thermodynamics for vacuum, stationary, asymptotically flat solutions with conical singularities 
reads 

dM = T H dS + ftdJ + PdA , M = M[S, J, A] . (2.1) 

The first two terms on the right hand side are the standard ones also present for regular solutions. 
These involve the Hawking temperature, Th , the entropy, S, the angular velocity of the horizon 
(or horizons), f2 and the angular momentum, J. In the presence of a conical singularity, which 
is exerting a pressure P, with world- volume spanning a spacetime area A/Th (computed in the 
Euclidean section), the last term should be included. Moreover, the energy functional M. is the 
ADM mass, Madm, subtracted by the energy associated with the conical singularity as seen by 
an asymptotic static observer, Ei nt . Generically, Ei nt = —PA Transforming to the energy 
functional Madm, corresponds therefore to a Legendre transform 

M ADM =M-PA, (2.2) 

leading to the first law in the form: 

dM ADM = T H dS + QdJ - AdP , M ADM = Madm [S, J, P] . (2.3) 

Observe that for the energy functional M. the independent variables are the extensive ones. 
That A, rather than P, is the extensive variable can be seen from an example. Consider the Z2 
invariant Israel-Khan solution with three Schwarzschild black holes. The conical singularity (and 
hence P) is the same in the two conical sub-sets, and the total A is twice that of each conical 
sub-set. It follows that thermodynamical equilibrium requires that Th,Q are the same in all 
connected components of the event horizon and that P is the same in all sub-spaces where there 
exist conical singularities. Of course, if there are no conical singularities, the two 'internal energy' 
functionals are equal, Madm = -M, and (|2.ip and (|2.3p coincide. But generically, M. plays a more 
fundamental role. This seems natural; indeed M, unlike Madm, takes into account the Komar 
masses of the individual black objects and the energy associated to the conical singularity. Thus, 
the remaining energy functionals (grand canonical potential, canonical potential and enthalpy) 
are obtained by the standard Legendre transforms from the energy functional AL 
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Alternatively, if the Legendre transform starts from the ADM mass, it must also include the 
variables associated to the conical singularity. For instance, performing the Legendre transform, 

W = M - T H S - QJ = M ADM - T H S — QJ+ PA , (2.4) 

we obtain the first law in terms of the grand canonical potential W (or Gibbs free energy) 

dW = -SdT H - Jdn + PdA , W = W[T H , 0, A] . (2.5) 

Thus, the entropy S, angular momentum J, pressure P, and energy functional M. of the system 
are given by 



S 



dW 



J = 



dW 



P = 



T H ,A 



dW 
~dA 



(2.6) 



In the Euclidean approach to black hole thermodynamics, the grand canonical potential is de- 
termined by 



W = T H I , 

where I is the Euclidean action. This action is 

Eint 



I = h- 



(2.7) 



(2.8) 



Iq is the contribution to the action from black objects found when neglecting the conical singu- 
larity (which arises from the boundary term). Consider the Euclidean action for the metric g 
over a region Y with a boundary dY to have the form Q : 



1 



16ttG 



R 



1 



d JY 



IttG 



d JdY 



(K - K ) , 



(2.9) 



where K (Kq) is the trace of the second fundamental form of the boundary dY (embedded in flat 
space) and Gd is Newton's constant in d— dimensions. 1 An asymptotically flat Euclidean metric 
with a single angular momentum parameter may be written in the form, near spatial infinity, 2 



ds z 



r d— 3 



dr l 



/' 



r d— 3 



dr 2 +r 2 {d6 2 +: 



cos 2 8dVL 2 d _ A ) 



2;/ sin 2 6> 



r d— 3 



where the parameters fj,, j relate to the ADM mass and angular momentum as 



16tt Madm 



" (d-2)n d _ 2 

£ld-2 denotes the area of the unit (d — 2)-sphere. 



Iil = 



8tt J a 



DM 



fid 



drdcj), (2.10) 



(2.11) 



1 In the following we take units in which Gd = 1. 

2 A higher dimensional black object may have k=[{d— l)/2] different angular momenta. Known examples with 
k > 1 are the d > 4 Myers-Perry black hole [j| and the d = 5 Pomeransky-Sen'kov double spinning black ring 
lldll . The generalisation of our results to higher k is straightforward. 
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The extrinsic curvature is denned as K^ v = — h^h^W a n^ , where h^ u = — n^n^. Taking 
dY to be the product of the time axis with a sphere of large radius tq , the unit normal to the 
hyper-surface dY is n p = dpT j \Jg rr . The traces we need are then 

K = gTK^ = d -l--J^ + o(^ l 

and Kq = (d — 2)/ro- Considering a space free of conical singularities; then the tree level 
Euclidean action will come entirely from the surface term, which gives 

Madm f , Madm a 

where (3 = 1 /Th is the periodicity of the Euclidean time. 
The second term in (|2.8|) is given by [l[ 

E mt = -PA = —A ; (2.12) 

it is the total energy associated to the strut as seen by a static observer placed at infinity. Here 
(5 is the conical deficit/excess associated to the conical singularity. Thus, the Gibbs free energy 
is always given by 

w(t h , a a) = M^M^hS _ i<? H n,A) A 

d — 2 07T 

It is now straightforward to compute the quantities (|2.6p . These will obey the formula 

d- 3 



_-M ADM =T H S + nj, (2.14) 

which is a generic consequence of relations (|2.2[) and (|2.4p . using also (|2.7p and (|2.8p . This is the 
Smarr formula; but note well that this is only the case because, in our description, S, J obtained 
from (|2.6p coincide with A#/4, Jadau where Ah is the sum of the areas of the event horizon 
connected components. 

We shall also compute, for the examples we consider, two other quantities, which are of rele- 
vance for the analysis of thermodynamical stability. The first one is the isothermal moment of 
inertia 



dJ 



(2.15) 

Th.A 



which is compute in the grand canonical ensemble. The second one is the specific heat at 
constant angular momentum, which is computed in the canonical ensemble. Thus we introduce 
the canonical potential (or Helmholtz free energy) via the Legendre transform 

F = W + QJ, (2.16) 

such that the first law takes the form 

dF = -SdT H + QdJ + PdA , F = F[T H , J, A] . (2.17) 
The specific heat at constant angular momentum is computed as 

dS 



Cj - TH df^ 



(2.18) 



J,A 
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The reason for working with Cj, rather than the specific heat at constant angular velocity 
Cfi = Tn{dS / dTn)\n ,a, is the following. In the grand canonical ensemble the condition for ther- 
modynamical stability is the positivity of the Weinhold metric, which amounts to the positivity 
of fiTff ,A an d Cj, rather than Cn (see @ and references therein). 

For completeness let us note that we can introduce a fourth thermodynamical potential, the 
enthalpy H, via the Legendre transform 

H = W + T H S. (2.19) 

The first law then reads 



dH = T H dS - JdQ + PdA , H = H[S, fi, A] . (2.20) 

We close this section by remarking that although the solutions of Einstein's equations with 
(naked) conical singularities should not be faced as vacuum solutions, this does not affect the 
generality of the relations derived above. Indeed, as in the simplest case of a cosmic string 
[ll| . the conical singularity is supported by a matter source with a precise form for its energy 
momentum tensor [l2fl . As discussed in [1[, however, the inclusion of the contribution of this 
matter source to the total action does not change the expression (|2.8p of the tree level action 
of the system. The reason is that taking into account the matter contribution, one should also 
subtract the contribution of a nontrivial background with an equivalent source, which consists 
in a finite piece of a string/strut (for d — 4) or a deficit /excess membrane (for d = 5) in a flat 
spacetime geometry. To see this explicitly, consider the total tree level action 

c ™)-i: (<f> K - i Ks) , (2.21) 

lay JdYs ) 



16tt J 8tt 



where the matter Lagrangian C m has been included and the reference background, wherein the 
trace of the extrinsic curvature is Kg, includes an equivalent source. The Lagrangian for the 
matter source that supports the conical singularity is C, n = —S/(8tt)Ss [U, where Ss is a Dirac 
delta function with support on the world- volume of the conical singularity. To compute the last 
term in (|2.2ip write the reference background as (for concreteness take d = 4) 

dsj = dr 2 + p 2 {SZ i oc + a 2 51 s + 5f)d<t> 2 + dp 2 + dz 2 , 

where 8 h a = 1 for z G [a, b] and zero otherwise. Then, considering dY to be the product of the 
times axis with a cylinder around the z axis, 

f Kg - f K a = 2tt/3 / (a - 1) dz = -5 Area, 

JdY 6 JdY J-z 

where S = 2tt(1 — a) is the conical deficit /excess present and Area = drdz is the area of the 
surface spanned by the conical singularity. Thus 

£ m + t- f K s = F" / Ko, (2-22) 

Y 87r JdY s 87r JdY 

which brings us back to the Euclidean action (|2.9[) . This argument reveals, again, the physical 
significance of the pressure P = —5/8ir and of its conjugate variable A = Area/f3; they define 
the physical properties of the matter source supporting the conical geometry. 

In the following sections we shall compute (|2.6[) . using (|2.13p . as well as (|2.15p and (|2.18p . for 
various examples. 
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3. BLACK RINGS ROTATING IN A SINGLE PLANE 



The rotating black ring solution in d = 5 Einstein gravity provides perhaps the simplest 
application of this formalism. The thermodynamics of the static solution was investigated in [l| 
where it was argued that, in the absence of rotation, all configurations are thermodynamically 
unstable. We shall now see that this conclusion still holds when rotation in either S 1 or S 2 is 
included. 



3.1. S 1 rotating black ring 

It was observed in [l3j that one can introduce rotation (along S 1 ) in the static black ring by 
performing a single soliton transformation; this is achieved by including a negative density rod 
on the rod structure of the seed metric to facilitate the addition of the angular momentum to 
the static black ring, when applying the inverse scattering method. The rod structures of the 
seed metric and of the resulting solution are shown in Fig. Q] The solution is characterised by 
four dimensionful parameters: the length of the two finite rods, i.e. a.32 and 043 (a^- = a; — a 3 ); 
the BZ (i.e. Belinsky-Zakharov, from the inverse scattering method) parameter 6; and the length 
<X2i of the phantom rod of the seed metric. However, it is useful to take out the overall scale of 
the solution and then introduce dimensionless parameters reflecting the length of the finite rods. 
We choose the overall scale L to be 

L = 04 — Q4 , 

and then introduce two dimensionless parameters t^i as 

Q-i — a l 

Ki = — — ; — , for i = l,2. 
L l 

Then we shifted the whole rod configuration along the z-axis, i.e. z — > z + ax, which explains 
the labelling of the rod endpoints in Fig. [T] (right panel). Due to the addition of the phantom 
rod, we are left with a singularity at (p = 0, z = 0). This shows up as a z" 1 divergence in the 
metric component g^ as indicated by the dots in Fig. [T](left panel). The singularity is, however, 
completely removed by setting 



b = ±V2^ 1 ~^L. (3.1) 

With this choice for b, the metric is completely smooth across (p = 0,2 = 0), and the final 
solution is fully characterised by L, Kx, and k 2 . The ex plic it line element can be found, in Weyl 
coordinates but in a different parameterisation, e.g. in [14| . 

The ADM mass, ADM angular momentum, angular velocity of the horizon O^,, Hawking 
temperature, Th and event horizon area, Ah, read: 

3tt 2 n 2 ^ 

Madm = —rL k 2 , Jadm = b , = — 



Tr _ 1 (l-«i) Ah _ - m)L 2 ( 3 - 2 ) 

2tt ^2k 2 (k 2 - Kx)L ' T H 

A straightforward computation shows that the physical quantities computed above satisfy the 
Smarr relation (|2.14p . if S = Ah /4 and if b is fixed according to (|3.1j) . 

Generically the solution contains a conical singularity. We choose it to be placed along the 
finite space-like rod. Observe that it may either be a conical excess S < 0, if the ring is under- 
spinning, or a conical deficit S > 0, if it is over-spinning. The excess/deficit is given by 

5 =4-^)- (m) 
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(i,o,o) (i,o,n) 



t i • 

(0,1,0) (0,1,0) (0,1,0) (0,1,0) 

«s 

(0,0,1) (0,0,1) ^ (0,0,1) 

a\ ci2 03 ^4 A4 x"2 1 



FIG. 1: Rod structure of the seed metric (left panel) and of the black ring solution, rotating along Si 
(right panel). The rods are located along the z-axis (z-axis) with p = in the left (right) panel. The solid 
rods have positive density and the dashed rod has negative density. The dots denote ring singularities 
at z = 0, which are removed by the fixing b according to (|3.1|l . 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

«2 «2 



FIG. 2: Parameter space of S 1 rotating black ring solutions. The regular solutions are along the solid 
line. On either side of this line they are not in mechanical equilibrium, being held by a membrane-like 
conical deficit (S > 0) or conical excess (8 < 0) (left panel). We also exhibit the sign of the mechanical 
moment of inertia of the solutions , together with lines of constant reduced spin (right panel). The dot 
marks the regular solution with 7 = 0; i.e. the meeting point between fat and thin ring branches. 



In Fig. [21 the sign of S is displayed, together with the mechanical moment of inertia, in 
parameter space. This mechanical moment of inertia reports the variation of angular velocity, fi 
with the angular momentum, Jadm, keeping fixed Madm and S, i.e. 



1 <9ft 



I OJadm 



1 — Kl 



7rL 4 ( — 4Kl + «2 + 3KlK2) 



The sign of / depends on ki and K2, but not on L. Thus, it can be completely exhibited by a 
plot on the ki — K2 plane - Fig[5J As usual, one can distinguish between the fat black ring branch 
with / > from thin black ring branch with / < 0, regardless of the value of 5. 

Requiring regularity on the rod < z < 1, fixes the period of <fi to be Atf> = 2tt, which 
corresponds to 8 = 0. From (|3.3p . one can see that this balance is achieved whenever 

K 2 = Kl(2 - Kl) . 



In particular, for black rings in equilibrium, / = if n\ = 2/3, which leads to the familiar reduce 
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0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 3: Sign of: specific heat at constant J, A, computed in the canonical ensemble (left panel); isother- 
mal moment of inertia, at constant A, computed in the grand canonical ensemble (right panel). Light 
(dark) grey regions correspond to positive (negative) sign. The solid line corresponds to regular black 
rings, wherein the ball separates the fat (below) from the thin (above) phases. Isotherms at constant A 
are also displayed in the right panel. 



spin value: 



•2 = 27^ J\ DM = Ki = 27 
■ 32 M\ DM k\ 32' 



where the thin and the fat balanced black rings meet. 

The parameter A, which is one of the relevant thermodynamical quantities, reads 



(l-^) 3/2 

- = an u 

P 1 - Kl 



Area „ r9 
A = — ^- = 2nL 2 



(3.4) 



where Area is the area of the surface spanned by the conical singularity, which is located at 
(p = 0, k,2 < z < 1). The grand canonical potential, from (|2.13p . is simply 



(3.5) 



It is not elegant to express in a simple way 3 the parameters L, K\, and K2 and thus W in terms 
of Th, and A. Nevertheless, when this last three quantities are used as thermodynamical 



3 For a balanced black ring, the Gibbs free energy has the simple expression Il5t 

1 



W 



i67r 2 r^ 
n 2 



(3.6) 
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variables, we precisely recover the Bekenstein-Hawking area law for the entropy: 



S = - 



dW 



H 



2 3 K 2 ~Kl r- -. r A H 

= 7T L°- V 2K 2(K2 - Kl) = — — . 

1 — Kl 4 



(3.7) 



Moreover, the thermodynamical angular momentum yields the ADM angular momentum and 
the pressure gives the expected relation (|2.12j) : 



J = — 



dW 

on 



Ja 



DM 



P = 



T H ,A 



dW 
dA 



T„,n 



5_ 

'Sir 



(3.8) 



Finally, observe that the thermodynamical mass A4, computed from Ai = W + ThS + 57 J, yields 
exactly (|2T2")l . upon using ([53]) - $£E§. 

The thermodynamical stability of the S 1 spinning black ring is analysed in Fig. [3J where the 
signs of the isothermal moment of inertia, 



ct h ,a = k 2 



4 (4-k 2 )k 2 + 4k 2 (2 + k 2 )-3kik 2 (4 + k 2 ) 



(4- k 2 )k 2 + Ki(-4 + k%) 
and specific heat at constant angular momentum, 



Cj = 7T L 



2 r 3 



iK 2 



(/t2 ~ ^i) 3/2 (^ 2 (8 + k 2 ) + «i(4 - 16k 2 + 3k 2 )) 
(ki - 1)((4 - k 2 )k 2 + 4k 2 (2 + k 2 ) - 3kik 2 (4 + k 2 )) 



are plotted. The plots are in the K\, k 2 space; indeed, as may be checked in the last two formulae, 
the L dependence does not change the sign of these quantities. It can be observed that there is 
no region in parameter space wherein both quantities are positive, and hence no region wherein 
the solution is thermodynamically stable in the grand canonical ensemble. This situation is 
analogous to that of the Kerr solution in four dimensions [l6[. In Fig. 01 isotherms for various 
values of A are plotted in the J-fl plane; the regions with negative and positive isothermal 
moment of inertia are clear and in correspondence with those plotted in Fig. [31 



3.2. S 2 rotating black ring 

This solution is explicitly given in [l7T | , where all relevant physical quantities where calculated. 
The metric is written in (t, x, y, <fi, ip) coordinates, where ip parameterises an S 1 and (a;, 4>) an S 2 . 
This particular black ring is rotating along the azimuthal direction </> of the S 2 . The coordinates 
x and y vary within the ranges —1 < x < 1, — oo < y < —1, with asymptotic infinity at x, y = — 1. 

After requiring regularity at infinity, the physical parameters for the black ring mass, angular 
momentum, temperature, angular velocity and horizon area are 



3ttR 2 X nR 3 \v Vl-A + ^ 2 (A- VA 2 -4i/ 2 ) 

Madm = -rr, r - ; — , J adm = — 7; r - — tTTTo > il = —- 



4(1 - A + v 2 ) ' _ (1 - A + 1/ 2 ) 3 / 2 ' 2XvR 
= v /(A 2 -4^)(^ 2 -l) ^ 2 B?\ 

H ^RX ' H (1 - A + v 2 )^/yj~l ' 

where 



(3.9) 



Uh 



-A + VA 2 - 4^ 2 
2^2 
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T «=- 




0.0 0.2 0.4 0.6 



FIG. 4: Isotherms - at constant Th,A - for the S 1 rotating black ring. Each (Jadm ,Th , A) generally 
corresponds to three values of Q, for low values of Th and A. For high values of Th and/or A, £t h ,a 
just changes sign once. 

determines the location of the horizon. The parameters A and v take values 2v < A < 1 + v 2 . In 
the limit v — > 0, we recover the static black ring. 

The solution contains a conical excess angle 6 along x = 1 and < y < — 1: 



S = 2n\l- 



1 + A + v 



1 - A 



(3.10) 



One can easily see that it is impossible to required reg ularity for nonzero A, which means that 
a S 2 rotating black rings cannot be in equilibrium |17| . Computing the mechanical moment of 
inertia 



9Jadm 



Mai 



(1 - A + v 2 f{\ 2 - Av 4 - AVA 2 - Av 2 ) 
27ri? 4 A 2 z/VA 2 - 4i/ 2 



it can be observed that just like for ordinary black holes, / is always positive - increasing Jadm, 
f2 also increases. 

The parameter A, related with the area of the surface spanned by the conical singularity, is: 



A 



Area 



ttR" 



1 + A + f 2- A- VA 2 



Av 2 



l-\ + v 2 + A + VA 2 - Au 2 



(3.11) 



The Gibbs free energy is given by (|2 . 13|) . From it, we precisely recover, using (|2.6[) . S = Ah/A, 
J = Jadm and P = —5/8ir. 

The thermodynamical stability of the S 2 rotating black ring is analysed in Fig. where 
the signs of the isothermal moment of inertia, €t h .a, and specific heat at constant angular 
momentum, Cj, are plotted. For any values of the parameters, when £t h ,a is positive, Cj is 
negative and vice- versa. Therefore, the S 2 rotating black ring is also thermodynamically unstable 
in the grand-canonical ensemble. 
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2.0 P 




FIG. 5: Sign of the isothermal moment of inertia, er H ,a, an d specific heat at constant angular momentum, 
Cj, in the parameters space of S 2 rotating black ring (left panel). Isotherms in the J-S7 plane, for various 
values of A (right panel). 

(l,fl) (1,X1) 
t 

(0,1) (0,1) (0,1) 

4 

ct\ a 2 «3 an 

FIG. 6: Rod structure of the asymptotically fiat, equal mass and angular momentum double-Kerr 
solution, obeying the axis condition. 

4. THE CO-ROTATING DOUBLE-KERR SOLUTION 

The approach proposed in this work can be extended to multi-black objects as well. For 
solutions with a non-connected event horizon, one needs to require that all connected components 
possess the same temperature and angular velocity, to ensure thermodynamical equilibrium. 
Interestingly, these conditions do not impose mechanical equilibrium and one can discuss the 
thermodynamics of multi-black holes with conical singularities. 

Perhaps the most obvious example here corresponds to the asymptotically flat d = 4 double 
Kerr solution, with two black holes having the same Komar mass and angular momentum. 
We further impose the axis condition 0, @, which guarantees that the ADM mass (angular 
momentum) is the sum of the Komar masses (angular momenta) of the individual black holes. 
The solution is then characterised by three physical quantities corresponding to the mass M and 
the angular momentum J of each black hole and the distance between them. 

The simplest way to parameterise the solution and the physical quantities is in terms of the two 
BZ parameters b and c, the length 021 = 043 — a of the rods associated to the black hole horizon 
and the length 032 = A of the rods between the black holes. The rod structure is represented in 
Fig. [SI Additionally, it is necessary to consider the following constraint @ 

(1 + 2c 2 + 6V) (-£_) 2 + 2c 1 ; b2c2 I = (1 - bcf , (4.1) 
\A + a J — c A + a 
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which guarantees that the axis condition is obeyed. The line element of the resulting configuration 
can be constructed from the data presented in 0, @| . 
The physical quantities are given by: 



Mi = M 2 

Jl = J 2 



(1 + b 2 ) (1 - be) (A + a) (X + 2a) j] 
2((l-bc)X + 2a) A ' 

(b(bc - 1)A + ((1 + b 2 )c - 2b)a) 

'(1 - bc)((l - &c) 2 A(A + 2a) + (1 - b 2 )(l - c 2 )a 2 f 
A ((6 2 - l)cA(A + 2a) + b(c 2 - 1)(A 2 + 2A<r + 2a 2 )) 



aMAr] 



T, 



2a(l + 6 2 )(A 2 +3Act + 2ct 2 ) 
A 2 



H 



A 1 - A 2 - — 
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ATtarf(\ + b 2 )(X 2 + '3Xa + 2a 2 ) 
a 

Th 



(4.2) 

(4.3) 

(4.4) 

(4.5) 
(4.6) 



where the individual masses and angular momenta are computed as Komar integrals and 
A= ((l-6c)A+(l + 6)(l-c)<r)((l-6c)A+(l-6)(l + c)<r) , v = {l-bc)\ + {l-c 2 )a . (4.7) 



It is now easy to check that once again the physical quantities satisfy the Smarr relation (|2.14p , 
taking Madm = M% + M2, the total angular momentum to be Jadm = Jl + J2 and the entropy 
to be (A$j + A%)/4. 

As argued by several authors 0, the spin-spin repulsion cannot balance the gravi- 

tational attraction in a double Kerr system where both Kerr objects are black holes. In this 
particular co-rotating limit one has the following excess angle along the section in between the 
black holes: 



6 

2^ 



(1 + Abc-c 2 -b 2 (l-c 2 ))a 2 



(1 - 6c) 2 A 2 + 2(1 - bc) 2 Xa - 4bca 2 
The parameter A associated to this conical singularity is 

A 



A 



A 



((1 - bc)X + 2a)((l - bc)X - 2bca) 



(4.8) 



(4.9) 



Although simple to calculate, the free energy of this co-rotating system, given by (|2. 13|) is 
already a quite long expression. A careful analysis, where the constraint (|4.1[) is cautiously 
considered, shows (numerically) the expected results. 



S 



dW 



H 



Ah 



A 2 



J 



dW 



J1 + J2 



p 



T H ,A 



dW 
~dA 



_5_ 

Sir 



(4.10) 

Note that this perfect match between the entropy and the sum of horizon areas, also between 
the different angular momenta, is not possible if instead of (|2.5p we consider that the variation 
of the free energy W is given by dW = —ScLTh — Jdfl — FdX as in @. 

The thermodynamical stability of the co-rotating double Kerr system is analysed in Fig. [3 
where the signs of the isothermal moment of inertia and specific heat at constant angular mo- 
mentum are plotted, for a particular window of the parameter space leading to physical solutions, 
i.e. with positive masses and horizon areas. Since Cj is proportional to A 2 , £t h .a is proportional 
to A 3 and the length A is always positive, the sign of the previous function just depends on (b, c), 
which explains the parameters space in Fig. [7] Again, as for the single Kerr and the single 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 

b 



FIG. 7: Sign of the isothermal moment of inertia, ct h ,a, and specific heat, Cj, in a window of the (6, c) 
parameters space for the co-rotating double-Kerr solution in thermodynamical equilibrium. Some curves 
of fixed (Th,A) are also displayed. 



rotating rings above, it can be observed that there is no region in parameter space wherein both 
quantities are positive, and hence no region wherein the solution is thermodynamically stable 
in the grand canonical ensemble. The isotherm behaviour, namely the variation of the angular 
moment with (ft, A) for fixed Tjj is in complete agreement with what is described in [6j. In 
particular the curve J = J(f2) is just like the one describing a single Kerr back hole, when the 
two black holes are far away or infinitesimally close. 



5. THE BLACK SATURN SOLUTION 

The black Saturn [l3[ is a d = 5 asymptotically flat solution describing a black ring around a 
concentric Myers-Perry black hole. Both objects have angular momentum only in a single plane 
- on the plane of the ring along the S* 1 direction. The solution was generated and thoroughly 
analysed in fl3j | , requiring regularity (on and outside the even horizons) . 

We are now interested in the particular limit where the solution is in thermodynamical equi- 
librium, i.e. when the black ring and the black hole have the same temperature and the same 
angular velocity, but not necessarily in mechanical equilibrium. The physical quantities of the 
general solution are explicitly given in [l3T | in terms of five parameters: L, Ki,rt2,/C3 and C2, 
where the last is a dimensionless parameter related to the original BZ parameter, L is the overall 
scale and the three dimensionless parameters m are related with the rod endpoints - Fig. [5] 
The physical quantities are explicitly presented in terms of this parameterisation in sec. 3.5 
and 3.7 of [13]. In what follows, we denote the black hole (black ring) mass, angular momen- 
tum, horizon angular velocity, temperature and horizon area by M BH , J BH ,ft BH ,T^ H , A^ H 
(M BR , J BR , n BR , T BR , A BR ). Observe that M ADM = M BH + M BR and J ADM = J BH + J BR - 

Setting VL BH = Q, BR and solving for C2 implies 



K 3 1 - Ki 

C2 = "j 

KlK 2 1 — K3 



(5.1) 
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(i,o,n BR ) (i,o,n BH ) 



(0,1,0) 



(0,1,0) 



(0,0,1) 



Af 3 K 2 K\ 1 



FIG. 8: Rod structure of the black Saturn solution. 



Considering the previous constraint and demanding T£ = leads us to define K3 in terms 
of K± and K2 through the relation: 

1 — ^2 

«3 = K\ + K 2 z ■ (5.2) 

1 — K\ — K 2 

It is important to note that the k, 's satisfy the ordering < K3 < n% < K\ < 1, and we can also 
add that (at least) for all physical solutions in thermodynamical equilibrium, K\ + k 2 > 1. 

The original solution has a conical singularity membrane in the plane of the ring, extending 
from the inner S 1 radius of the black ring to the horizon of the S 3 black hole. Applying the 
constraints (|5.ip and (|5.2p that lead to equilibrium thermodynamics, the corresponding deficit 
angle is 

8_ k 2 (1 - k 2 ) 2 I K2 

2lT («i - «2)(«1 + «2 — 1) V «l(«l«2 - («1 + k 2 - l) 2 ) 

where we can see that it is still possible to reach mechanical equilibrium, taking into account the 
properties of k^'s. In Fig. [H]the sign of 8 is analysed along the parameter space of black Saturn 
solutions in thermodynamical equilibrium. The balanced solutions lie along the solid line and 
have already been study in [2l| . This solutions can also be classify as fat or thin in function of 
the sign of the mechanical moment of inertia, /, just like single black rings. 

The sign of / is displayed in right panel of Fig. [SJ just in function of (k\,K2) because / is 
proportional to L 4 and so L does not influence its sign. Analysing, along the parameter space, 
the family of solutions with constant reduced spin 

, 2 _27tt (J br + J bh ) 2 
Jbs _ 32 {M BR + M BH f ' 

one can observed that, as expected, 8 < for under-spinning back Saturn configurations, while 
for over-spinning ones, 8 > 0. Generically, most of the angular momentum is in the ring. In the 
fat branch, which is associate with lower reduce spin values, the fraction of the total mass that 
goes into the central black hole increases with k 2 . For high jss, the ring also carries most of the 
total mass. 

The parameter A is 

A = 2-kL . — (/ci« 2 - («i + K 2 - I) 2 ) , (5.3) 

K 2 (l - K2) Z V K 2 

and the Gibbs free energy of this system is given by 
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FIG. 9: The coloured region represents the space of black Saturn solutions in thermal equilibrium. 
The regular solutions are along the solid line. On either side of this line they are not in mechanical 
equilibrium, being held by a membrane- like conical deficit (5 > 0) or conical excess (5 < 0) (left panel). 
We also exhibit the sign of the mechanical moment of inertia of the solutions and some lines of constant 
reduced spin (right panel). 




FIG. 10: Analysis of the black Saturn solution in thermodynamical equilibrium for fixed ADM mass and 
reduce spin jss = 1- The plots exhibit how the relative mass M BH /M BR varies with the reduce spin of 
the black ring jbr = \/27ir /32J BR /M^£i M (left panel), the reduce angular velocity fi = \J 8A1ad m / 3ir£l 
(middle panel) and the reduce total horizon area Ah = \/27/(256ttMj, dai )(A h r + A H H ) (right panel). 
The arrow indicates the variation of these quantities along the line j B g = 1 in Fig[9] from K2 = to 

Kl = 1. 



taking into account the constraints (|5.1[) and (|5.2[) . 
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The appropriate set of thermodynamical variables is (Th , ^2, A), where 



T H = 



2-kLkx 



1 



1 — K\ (ki + k 2 — l)(«i — k\ + n 2 — K 2 ) 

1 



(5.5) 



^/(2(«i + k 2 - 1)(kiK2 - (ki + «2 - l) 2 )) ' 
is the temperature of the system that results from Th = = T^ H ; 

(1 - K 2 )K 2 / K 2 (kiK 2 - «1 + ti\ - K 2 + K 2 ) 



L(«l + K 2 — 1) V 2/Cl(«l + K 2 — 1)(«1 — K 2 + K2 — «1 



2\2 : 



(5.6) 



is angular velocity of the black ring and the black hole, i.e. ft = £l BR = il BH , and .4 is defined 
in (|5.3p . The entropy is 



ST, 



(I _ K . )2 

tt 2 L 3 )- (5-3 («i + k 2 - l)(2«i + k 2 - l)(«i - k 2 + k 2 - k|) 

(1 — K2) K 2 



(1 - k 2 )(1 - 1k,\ + K'f - 2k 2 + K1K2 + k\) 1 



1(A BR + A BH ) ■ 

"V (l-«i)(«i + «a-l) "4 ( ^ ); 

so, once again, the Bekenstein-Hawking area law is recovered. A similar calculation shows that 



J = 



dW 



= J BR +J BH 



p 



T H ,A 



dW 
ThA 



_6_ 

'8tt' 



From the above results, one may again conclude that the thermodynamical mass M, computed 
from M = W + T H S + OJ, yields exactly (|Z2f 



M = M BH + M BR -^-A = M ADM 



Finally, the thermodynamical stability of the black Saturn is analysed in Fig. [TTJ where the 
signs of the isothermal moment of inertia and specific heat at constant angular momentum are 
plotted. Again, as for the single and double Kerr and the S 1 rotating ring, there is no region in 
parameter space wherein both quantities are positive, and hence no region wherein the solution 
is thermodynamically stable in the grand canonical ensemble. 



FINAL REMARKS 



In this paper, we have further developed the proposal in [l| for the thermodynamical descrip- 
tion of asymptotically flat solutions with conical singularities. The main observation in [l| is 
that this description yields the Bekenstein-Hawking formula for the entropy which is obtained 
by differentiating the Gibbs free energy. In the Euclidean approach to quantum gravity, the 
Gibbs free energy is obtained from the Euclidean gravitational action. Observe that in previous 
approaches, deviations from the Bekenstein-Hawking formula were obtained, in solutions with 
non-connected event horizons (22|. Here, we have considered various examples of stationary so- 
lutions, and showed that the description gives the natural results; thus, the thermodynamical 
angular momentum, obtained again by differentiating the Gibbs free energy coincides with the 



17 




FIG. 11: Left panel: sign of the specific heat. Right panel: sign of the isothermal moment of inertia. 



ADM angular momentum, again in contrast with previous descriptions Q. Given these results, 
which support the idea that we have a reliable thermodynamical description of these solutions, 
we have also considered their thermodynamical stability. We have found that there is no point 
in the parameter space for which the appropriate Hessian matrix is positive definite. Thus these 
solutions are always unstable in the grand canonical ensemble. 

It is worthwhile remarking that, similarly to the static case in [l|, the location of the conical 
singularity was a matter of choice. After a suitable re-scaling, all solutions considered in this 
work have an alternative interpretation as non-asymptotically flat black objects (i.e the coni- 
cal singularity may be chosen to extend to spatial infinity). The asymptotic spacetime then 
corresponds to a cosmic string spacetime for d = 4, or to its higher dimensional analogue (a 
membrane for d = 5). The conical deficit /excess 5°° differs from 6 of the corresponding solution 
with the conical singularity having a compact support in the bulk, the relation between these 
two quantities being 




The formalism proposed in Section [5] can easily be generalised to this situation and the ther- 
modynamical behaviour of the system, in particular its instability, should be independent of the 
choice for the location of the conical singularities. 
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FIG. 12: Isotherms - at constant Th,A - for the S 1 rotating black ring. Each (J, Th,A) generally 
corresponds to three values of Q for low values of Th and A. For high values of Th and/or A, £t h ,a 
just changes sign once. 
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